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Understanding double-resonant Raman scattering in chiral carbon nanotubes: 
Diameter and energy dependence of the D mode 
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We present a theoretical model to describe the double-resonant scattering process in arbitrary 
carbon nanotubes. We use this approach to investigate the defect-induced D mode in CNTs and 
unravel the dependence of the D-mode frequency on the CNT diameter and on the energy of the 
resonant optical transition. Our approach is based on the symmetry of the hexagonal lattice and 
geometric considerations, hence the method is independent of the exact model that is chosen to 
describe the electronic band structure or the phonon dispersion. We finally clarify the diameter 
dependence of this Raman mode that was controversely discussed in the past and demonstrate that, 
depending on the experimental conditions, in general two different dependencies can be measured. 

We also prove that carbon nanotubes with arbitrary chiral index can exhibit a D mode in their 
Raman spectrum, in contrast to previous symmetry-based arguments. Furthermore, we give a 
direct quantification of the curvature-induced phonon frequency corrections of the D-mode in carbon 
nanotubes with respect to graphite. 

PACS numbers: 61.48.De,78.67.Ch,63.22.Gh,78.30.Na 


I. INTRODUCTION 

Carbon materials such as graphene and single-walled 
carbon nanotubes (CNTs) attracted much scientific in¬ 
terest in the past decades^’^. Their extraordinary prop¬ 
erties made them a perfect candidate for applications in 
novel electronic devices such as transistors or sensors^^^. 
Here, functionalized systems are most important, since 
their properties can be designed and manipulated as 
preferred®^®. The success of a functionalization is of¬ 
ten monitored by investigating the increase of the Ra¬ 
man D-mode intensity®. This mode stems from TO 
phonons around the K points of the graphene Brillouin 
zone and is activated by breaking the translational in¬ 
variance, thus relieving selection rules^®. Although the 
defect-related origin of the D mode was known for a 
long time, the underlying Raman processes remained 
unclear^^. Thomsen and Reich solved this longstand¬ 
ing question by introducing the double-resonant scatter¬ 
ing theory for graphitic systems^^, explaining the anoma¬ 
lous dispersion of certain Raman modes, e.g., the D and 
2D modes, with laser excitation energy^^. However, the 
double-resonant D mode was thought to exist only in car¬ 
bon nanotubes that satisfy kp = where fci? is the 
reciprocal vector of the singularity in the joint density 
of states. This relation is particularly true for armchair 
carbon nanotubes and other metallic tubes with TZ = Z, 
i.e., (ni — n2)/(3n) = integer^’^®. In contrast, semicon¬ 
ducting carbon nanotubes were predicted to not exhibit 
a D mode with the same systematic excitation energy 
dependence of its frequency^^. However, recent exper¬ 
iments on semiconducting CNT samples enriched with 
single chiral indices demonstrated the existence of a D 
mode in CNTs that do not satisfy the aforementioned 
restriction^®. Furthermore, there is an ongoing discussion 
about the dependence of the D-mode frequency in CNTs 
on the tube diameter and the transition energy. During 


the past decade many different explanations were pos¬ 
tulated and controversely discussed. For instance, Souza 
Filho et al}"^ proposed an inverse diameter dependence of 
the D-mode frequencies, whereas Refs. 18 and 19 claimed 
a dependence proportional to the tube diameter. Thus, 
the systematics of these scattering process is still not fully 
understood and needs clarification. 

In this work, we derive a universal, geometrical model 
that describes the dependence of the resonant phonon 
wave vector in the double-resonant scattering process on 
the tube diameter and the transition energy of the pro¬ 
cess for all tubes with arbitrary chiral indices, rejecting 
the previously predicted restriction of the D mode to par¬ 
ticular CNTs. Our model is based on the symmetry of 
the hexagonal lattice only and is thus universally valid, 
no matter which approximation for the electronic band 
structure or phonon dispersion is used. We apply this 
model to investigate the diameter dependence of the D 
mode in CNTs for the resonant optical transition, which 
is still under controverse discussion. We will point out 
that, depending on the specific experimental conditions, 
in general two different diameter dependencies can be ob¬ 
served. Furthermore, we will highlight the importance of 
nanotube curvature effects on the phonon frequencies to 
explain the experimentally observed dispersion of the D 
mode. Finally, we derive a diameter correction for the 
D-mode frequencies in carbon nanotubes with respect to 
the D-mode frequency in graphite. 


II. SIMULATION DETAILS 

All geometric considerations in this paper are based 
on the symmetry operations of the hexagonal lattice of 
graphene. Hence, within the framework of zone-folding, 
the results do not depend on the choice of the model 
to describe the electronic band structure or the phonon 
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FIG. 1. Illustration of the double-resonant scattering process 
between two energetically equivalent minima in the helical 
band structure of a (5, 3)-tube. The transition energy is la¬ 
beled as Eii, the position of the electronic minimum as fcmin, 
and the resonant phonon wave vector as q^es- The transla¬ 
tional period is denoted with a; q reflects the subbands in 
linear quantum numbers. 



FIG. 2. Illustration of the dependence between the electronic 
transition and the resonant phonon wave-vector for a coprime 
tube, i.e., n = 1. The black solid line denotes the helical vec¬ 
tor of a (3, 2)-tube, the full orange circles mark a minimum 
in the electronic band structure (femin) and the correspond¬ 
ing resonant phonon wave-vector length (^res), which is twice 
fcmin. The open red circles denote the closest K points to 
fcmin and qres, respectively. As can be seen easily, the dis¬ 
tance h — K\ — fcmin is always twice the distance K 2 — qies- 


dispersion. 

The calculcations of electronic band structures and 
phonon dispersions in this work are based on the POL- 
Sym code in the sixth-nearest-neighbor approximation^*^. 
This package uses the modihed group projector technique 
and includes curvature effects in the calculations of the 
band energies and phonon frequencies. Electronic bands 
and phonon dispersions are calculated for all 274 chiral 
tubes with diameters between 5 and 25 A. Furthermore, 
we used an experimental graphite phonon dispersion from 
Ref. 21 for an alternative calculation of the H-mode fre¬ 
quencies and for comparison with the POLSym-derived 
frequencies. 

We model the double-resonant scattering process in 
carbon nanotubes by assuming that scattering occurs 
mostly between equivalent minima in the electronic band- 
structure (compare Fig. 1)^®. As demonstrated by previ¬ 
ous works, the double-resonance process is dominated by 
the incoming resonance with excitonic transitions at van- 
Hove singularities in the electronic bandstructure^^’*^*^, 
justifying our assumption. Furthermore, we assume that 
the phonon frequency /iWphonon = OmeV, as the reso¬ 
nant phonon wave-vectors in the double-resonance pro¬ 
cess depend only weakly on the energy of the involved 
phonons^'^. These assumptions do not affect the general 
validity of our results in this work. 


III. RESULTS AND DISCUSSION 

As we mainly consider chiral carbon nanotubes, it 
is convenient to use the helical representation of the 
wave vector and electronic bands for the following 
analysis^^’^®. In the zone-folding approach, the electronic 
band structure in helical quantum numbers is given by 
the cut of the helical wave vector with the electronic 
TT and TT* bands of graphene along its path in reciprocal 


space. The helical wave vector is given by 

= - (-712 ki -I- ni k2) + Aki(ni,n2) + Ak2(ni,n2), 
n 

where ni and 712 are the chiral indices of the CNT and n 
is the greatest common divisor of (ni, 712 ). ki and k 2 are 
the unit vectors in graphene’s reciprocal space (compare 
Fig. 2). The Ak^ describe the distance in reciprocal space 
of the fth subband from the subband containing the F 
point of graphene and is non-zero only if ti > 1. For a 
more detailed derivation of the helical wave vector kj, and 
an introduction to the concept of helical zone-folding we 
refer the reader to the Supplementary Material^^. 

The D mode results from a double-resonant intervalley 
scattering process, including a transverse optical (TO) 
phonon and a defect *^^44^ jjj Figure 2 we show the sys- 
tematics of the D-mode scattering process for a (3, 2)- 
tube. Our approach is of course valid for all other tubes 
having more than one subband (see Supplementary Ma¬ 
terial). Without loss of generality, we assume that the 
minimum in the band structure of a CNT occurs at the 
position where the helical vector crosses a. K — M — K' 
high-symmetry line. The minimum shall have the k- 
vector fcmin (measured from the F point), hence the reso¬ 
nant phonon wave vector in the D-mode scattering pro¬ 
cess has the length This implies the so- 

called (7res ~ 2 Amin rule. We label the K point closest to 
the electronic minimum as Ki and refer to its distance to 
kmin as h. In general, the distance h is a measure for the 
transition energy, ie., a smaller value for h means a lower 
energy and vice versa. The three points F, Amin, and Ki 
form a triangle, which is indicated by the dark-gray area 
in Fig. 2. By similarity transformation of this triangle, 
we obtain the light-gray triangle formed by the points F, 
Qres) and 7^2. It is an intrinsic property of the hexagonal 
lattice that twice the distance between F and a K point, 
is again a K point ( 7 ^ 2 ). If Ki is the closest K point 
to Amin; then K 2 is the closest K point to gj-es- Since 
all sides of the larger triangle were doubled, the distance 
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FIG. 3. (a) Distance between the resonant phonon wave-vector and the closest K point as a function of the tube diameter 
for the resonantly enhanced scattering process obtained by a sixth-nearest neighbor tight-binding model for electronic band 
structures and phonon dispersions, (b) Calculated D-mode frequency a;D as a function of the tube diameter, (c) Calculated 
D-mode frequency as a function of the laser energy. The experimentally observed jumps between the different transitions are 
indicated. 


K 2 — Qres is now twice the distance h. This means that 
an electronic transition with a distance h from Ki results 
in a resonant phonon vector with a distance 2h from K 2 
(at the same angle with respect to the closest K point). 
Thus, h and 2h are indicative for the energy of the tran¬ 
sition and the D-mode frequency, i.e., a larger value 2h 
in general means a larger TO phonon frequency and vice 
versa. From this we can deduce two important results: 
First, tubes with large diameters and thus lower transi¬ 
tion energies exhibit a systematically lower D-mode fre¬ 
quency for the resonantly enhanced transition compared 
to small-diameter tubes. Second, different transitions in 
a CNT exhibit different D-mode frequencies. Although 
these results seem trivial at a first glance, they enable a 
deeper understanding of the D mode in CNTs. 

If we now also consider trigonal warping effects in the 
electronic bandstructure, the electronic minimum may 
not be exactly on the K — M — K' high-symmetry di¬ 
rection, but slightly shifted away. This would lead to 
a slightly different distance h and angle between femin 
and Ki. Nevertheless, the distance between the resonant 
phonon wave vector q^es and K 2 is again 2h. Hence, this 
relation is independent from the exact position of the 
electronic minimum, which may differ depending on the 
model used in bandstructure calculations. Thus, our re¬ 
sults are universally valid as long as zone-folding is an ap¬ 
propriate model to describe the properties of CNTs^®“^°. 

Figure 3 (a) presents the calculated distance between 
the resonant phonon wave vector and the closest K point 
for different optical transitions as a function of the tube 
diameter d. As can be seen, the distance 2h decreases 
with a 1/d-dependence. Ergo, the resonantly enhanced 
TO phonon frequency decreases likewise. Furthermore, 
energetically higher transitions have a larger distance 2h 
for the same tube diameter, corresponding to higher D- 
mode frequencies. Since the distance 2h directly depends 
on the energy of the optical transition of a tube, we 
observe a close correspondence between the dependence 
shown in Fig. 3 (a) and the so-called Kataura plot^^. By 


translating the distance 2h into a phonon frequency us¬ 
ing the POLSym-calculated phonon dispersions, we ob¬ 
tain the D-mode frequencies shown in Fig. 3 (b). As dis¬ 
cussed before, a larger value for 2h reflects a larger TO 
phonon frequency and vice versa, thus, the systematics 
of this plot again resemble the Kataura plot. We observe 
a decreasing D-mode frequency with increasing tube di¬ 
ameter. Furthermore, electronic transitions with a higher 
energy show an increased frequency. 

If we plot the D-mode frequency as a function of the 
transition energy in Fig. 3(c), we observe an interest¬ 
ing feature in the D-mode dispersion. The dispersion of 
the branches is not continuous. In fact, a small jump 
in frequency between different electronic transitions is 
observed. This frequency jump has been reported pre¬ 
viously in resonance Raman studies on enriched carbon 
nanotube samples^^’^^’^^. As we discuss below, the ob¬ 
served jump is a direct consequence of curvature effects 
on the phonon frequencies and thus cannot be explained 
using a simple zone-folding approach. Considering a 
zone-folding third-nearest-neighbor tight-binding model 
without any curvature effects, only the helical vector 
would depend on ni and n 2 , but not the electronic bands 
nor the phonon dispersion. Hence, two different tubes 
with different diameters but the same transition energy, 
i.e., same distance h, would exhibit the same D-mode 
frequency, i.e., same distance 2h. Thus, by neglecting 
curvature effects, the experimentally observed frequency 
jump cannot be explained. However, it is well-known 
that phonons in carbon nanotubes show a strong depen¬ 
dence on the tube diameter due to curvature effects and 
rehybridization of cr and tt orbitals. For instance, in semi¬ 
conducting tubes, the F-point TO phonon frequency {G~ 
mode) decreases with decreasing tube diameter^^; a sim¬ 
ilar diameter dependence of the TO branch is expected 
around the K point. Thus, two different tubes with dif¬ 
ferent diameters but similar transition energies will not 
have the same D-mode frequency, i.e., the smaller tube 
has a lower frequency compared to the larger tube. Ergo, 
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FIG. 4. (a) Schematized experimental data of the D mode 
in carbon nanotubes (after Refs. 32 and 33). The solid lines 
denote fits to the experimental data for each transition; the 
numbers at each line correspond to the diameter of the tubes 
(in units of A). The dashed line indicates the D-mode dis¬ 
persion in graphite, (b) D-mode frequency correction as a 
funtion of the tube diameter. The red, dashed line reflects a 
fit to the data points. Solid symbols were obtained from (a), 
open symbols are single-tube measurements from Refs. 32 
and 33. (c) Curvature-corrected D-mode dispersion from ap¬ 
plying Eq. (1) on the graphite dispersion. The experimentally 
observed discontinuity could be reproduced. 


in a resonant measurement, one would observe a jump 
between different transition branches. In our calcula¬ 
tions, this discontinuity is smaller than experimentally 
observed. Furthermore, D mode from the Eii transition 
shows nearly no dispersion with transition energy. We 
attribute these discrepancies to the fact that the Kohn 
anomaly and the region around the K point in the TO 
phonon dispersion is not well approximated by the model 
for the phonon dispersion. Since the resonant phonons 
from the Efi transition stem from a region close to the K 


point, small deviations in the phonon dispersion directly 
influence the results and, in our case, lead to an over¬ 
estimated D-mode frequency with nearly no dispersion. 
We expect a larger jump and a better correspondence to 
the experiments for calculations using a non-orthogonal 
tight-binding model or an ab-initio approach. Neverthe¬ 
less, the existence of a jump can be reproduced in our 
simulations, as seen between Efi and E 22 in Fig. 4 (c). 

Next, we want to investigate the diameter distribu¬ 
tion along a given optical transition Ea . The energy of a 
transition is inversely proportional to the CNT diameter, 
hence, when following a D-mode branch while decreas¬ 
ing the transition energy [see Fig. 3(c)], the diameter 
decreases. 

The above presented results finally harmonize the dif¬ 
ferent conclusions from various previous works where a 
decreasing D-mode frequency for a decreasing tube diam¬ 
eter was claimed^®’^® or vice versa^^. It is very important 
to point out that, depending on the experimental condi¬ 
tions, one can generally observe two different dispersion 
behaviors. By measuring the same optical transition for 
all tubes, e.g., the E 22 transition, an cxl/d-dependence 
will be observed due to the mentioned diameter distri¬ 
bution along the E 22 branch. This would correspond to 
an experimental setup where the laser energy is tuned 
to measure every tube always at its resonance energy. 
By using only a single or just very few laser energies, a 
different behaviour is observed^®’This experimental 
condition corresponds to a vertical cut through the dis¬ 
persion relation shown in Fig. 3 (c). Along such a cut, a 
higher D-mode frequency corresponds to a higher tube 
diameter and a dependence proportional to d is observed. 
Both cases are distinctly different and must be separated 
carefully. 

To finalize our results, we want to give an explicit ex¬ 
pression for the correction of the D-mode frequency as 
a function of the tube diameter, as a correction to the 
D mode in graphite. As pointed out before, the phonon 
frequencies in CNTs show a strong dependence on the 
tube diameter. This dependence was shown both exper¬ 
imentally and theoretically for the as well as 

for the G“ mode®"*. However, a direct evaluation of this 
dependence for the D mode is still lacking. Starting from 
the schematized experimental D-mode dispersion shown 
in Fig. 4 (a) taken from Refs. 32 and 33 and the exper¬ 
imentally known dispersion of the D mode in graphene 
and graphite^°’®®, we calculate the frequency difference 
between both dispersions and plot them as a function of 
the tube diameter. Since the D mode in graphene and 
graphite does not include any curvature effects, the differ¬ 
ence between those dispersions is a direct measurement 
of the diameter-dependent correction of the D-mode fre¬ 
quency. The result of this calculation is shown by filled 
and open circles in Fig. 4 (b). As expected, the frequency 
correction is always negative and converges to zero for in¬ 
creasing tube diameter. The dashed line reflects a fit to 
the data derived from a classical mechanical model (see 
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Supplementary Material), using the formula: 


AuJoid) = A 





- 1 


quence of the diameter dependence of phonon frequencies 
in carbon nanotubes. 


( 1 ) 

IV. CONCLUSION 


Here, d is the tube diameter and H is a fit parameter 
obtained from fitting the data points in Fig. 4 (b) and 
was determined as H = 593 cm“ ^. 

In the following, we will apply the above result to cal¬ 
culated I?-mode frequencies obtained from a phonon dis¬ 
persion that does not include curvature-induced effects. 
We use the calculated electronic minima and resonant 
phonon wave vectors from our sixth-nearest-neighbor 
model and obtain the resonant phonon frequencies from 
an experimentally measured TO phonon dispersion of 
graphite from Ref. 21. The CNT diameters in this cal¬ 
culation were chosen to fit the experimentally measured 
nanotubes from Refs. 32 and 33, i.e., only tubes with di¬ 
ameters between 6 A and 14 A were considered. Without 
curvature corrections, a linear dependence between tran¬ 
sition energy and U-mode frequency is observed. By ap¬ 
plying the diameter-dependent frequency correction ac¬ 
cording to Eq. (1), CNTs with smaller diameter observe 
a larger frequency downshift than nanotubes with large 
diameters with respect to the graphite U-mode disper¬ 
sion. Thus, the slope of the D mode dispersion in CNTs 
is reduced compared to the dispersion in graphite [com¬ 
pare Fig. 4 (a)]. As a consequence, a discontinuity opens 
between different transitions. In our curvature-corrected 
U-mode dispersion from Fig. 4 (c), this discontinuity can 
be seen again most clearest between and E 22 (ap¬ 
prox. 15cm“^). For energetically higher transitions, the 
frequency jump decreases, e.g., the difference between 
EI 2 and Ef( is approx. 9 cm Although the calcu¬ 
lated frequency difference is smaller than experimentally 
observed, the experimentally observed discontinuity is re¬ 
produced in our simulations and is shown to be a conse¬ 


In summary, we derived a theroretical model to de¬ 
scribe the double-resonant scattering process in arbitrary 
CNTs. We applied this model to describe the diameter 
dependence of the D mode in carbon nanotubes for the 
resonant optical transitions. The presented approach is 
independent of the specific calculational model used for 
the electronic band structure or phonon dispersion and 
is therefore universally valid. We showed that, depend¬ 
ing on the experimental conditions, in general two dif¬ 
ferent diameter dependences can be observed. Further¬ 
more, we proved that the experimentally observed dis- 
contuinity between different transition branches is due 
to curvature effects that alter the TO phonon disper¬ 
sion at the K point. Finally, we derived a quantifica¬ 
tion of the diameter-dependent frequency correction of 
the D mode in carbon nanotubes with respect to the D 
mode in graphite. The presented results answer the long¬ 
standing question regarding the diameter dependence of 
the D in carbon nanotubes. Finally, we want to explic¬ 
itly point out that our approach is equally valid for all 
other double-resonant Raman modes in CNTs and can 
be easily expanded to intravalley scattering processes. 
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Helical wave-vector backfolding 


Since this paper mainly considers chiral carbon nanotubes, it is more convenient to use the helical representation of 
the wavevector and electronic bands instead of the conventional linear representation^’^. Therefore, we first give a 
short introduction and motivation to the helical quantum numbers. We will show that the electronic band structure 
in helical space is given by the cut of the helical wavevector with the electronic tt and n* bands along its path in 
reciprocal space. 

In linear quantum numbers, the electronic bands are given along the quasi-continuous reciprocal vector k^, which has 
the form 


k. = 



( 1 ) 


Here, ki and k 2 are the unit vectors in reciprocal space [compare Fig. 1 (a)] and ni and n 2 are the chiral indices of 
the tube. The number of graphene hexagons per unit cell is q, which also gives the number of electronic subbands^. 
In linear space the vector k^ is restricted to the interval (— 7r/a, tt/o], where the brackets ( and ] indicate open and 
closed interval limits, respectively. Here, a is the translational period along the tube axis. The electronic dispersion 
of the q subbands is then given by q cuts of k^ with the tt/tt* bands of graphene. However, q increases rapidly even 
for non-high-index tubes, e.g., q = 386 for a (9,7) tube. Therefore, the linear representation becomes complicated 
very quickly. The transformation from the linear to the helical representation resolves this problem and is therefore 
much more convenient for displaying band structures or phonon dispersion relations^. 

If kj, is transformed into its helical representation k^,, the interval is extended by the factor q/n, where n is greatest 
common divisor of (ni,n 2 ).^ This means that the q subbands in linear space are “rearranged” to form n subbands in 
helical space. Thus, the helical vector is defined in the interval {—qw/na, qn/na] and has the form^ 

k„ = - {-712 ki -I- ni k2) (2) 

n 

As before, the electronic band structure of CNTs in the helical representation is given by the n cuts of the helical 
vector with the tt/tt* bands. The different subbands are labelled with the helical quantum number m. In contrast to 
the linear represenation, the prefactors of ki and k 2 in the helical image are by defintion always integers, simply ni 
and n 2 - This means that the m = 0 subband always connects two F points on the reciprocal hexagonal grid, which are 
given by F = (0,0) and Fi = (—n 2 /n, ni/n) in units of the reciprocal lattice vectors. Hence, the length of the helical 
vector depends solely on ( 711 , 712 ) and is therefore directly proportional to the diameter of a carbon nanotube. This 
result is depicted in Fig. 1 (a) for two different chiral tubes. It is more convenient to use the helical representation 
of the electronic bands, especially for coprime tubes (n = 1), since all information about the electronic dispersion is 
given in just one subband. For tubes with more than one electronic subband {n > 1), the other subbands can be 
found by performing a parallel translation of the helical vector. A detailed analysis how of to obtain the shifted m 7 ^ 0 
subbands can be found in Ref. 1. 

Next, we want to introduce the concept of backfolding the helical band structure of carbon nanotubes. This approach 
is closely related to the usual sampliing og the Brillouin zone that is already known for linear bands in CNTs^’®, but 
was never applied to their helical bands. In Fig. 1 (a) we illustrate the folding of helical bands for two CNTs: for a 
(5,3)-tube, where k^ runs close to the armchair direction - thus introduced here as an “armchair-like” tube (ACL) - 
and for a (4, l)-tube, named “zig-zag-like” (ZZL) accordingly. The helical band of an ACL tube can now be folded 
back along the dashed lines running through F points in zig-zag direction, and vice versa for ZZL tubes along those in 
armchair direction. Using this approach, the helical vector along F — Fi (black line inside the blue shaded rectangle) 
can be reduced to the sawtooth-like path along F — r 2 (red line inside the red shaded rectangle). While both paths 
are of course equivalent, it becomes apparent that the amount of hexagons covered by the helical vector reduces 
significantly and the location of minima becomes even more obvious. 
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Figure 1. (a) Backfolding of the helical band structure of single-walled carbon nanotubes. The band structure of the “armchair- 
like” (5, 3)-tube is projected onto the zig-zag direction, whereas the “zig-zag-like” (4, l)-tube is projected onto the armchair 
direction. The band structures along the F-Fi-line can be represented by the folded path along F — r 2 . For the (4,1) tube, we 
also show the linear bands (green), (b) Helical band structure of (99, 97)-tube (gray, background) and a (5, 3)-tube (thick lines, 
foreground), q is the number of hexagons per unit cell, n Is the greatest common divisor of (ni,n 2 ), and a is the translational 
period of the tube along the axis. 


The consequences of the backfolding become more clear, when we consider the limit of large-diameter tubes, e.g., 
a (99,97)-tube. As shown before, the length of the helical vector is directly proportional to the tube diameter. 
Consequently, the helical vector gets folded at an increased number of positions for large tubes. In the following, we 
want to refer to this folding as the sampling of the Brillouin zone. Since for the ACL tubes (5, 3) and (99, 97) n 2 equals 
ni — 2, in both cases the zig-zag component of their helical vector will span cross two “rows” of hexagons and thus 
the projection area for both tubes is the same (see Fig. 1). Hence, the folded path in the red-shaded rectangle must 
cover the same three underlying hexagons more densely for the (99, 97)-tube compared to the (5, 3)-tube. Naturally, 
this means that the electronic tt bands of the graphene analogon are simply sampled more closely. This effect can be 
nicely seen in Fig. 1 (b), which shows the helical band structure of both tubes along their normalized helical vector. 
Here it becomes obvious that the band structure of the (99,97)-tube simply corresponds to a denser sampling of the 
same area as in case of the (5, 3)-tube. Since both projection areas are equivalent, both band structures could also 
be imagined encased in an evenlope function representing the band structure of the three graphene hexagons in their 
projection area. In our example, there are three F, M and K points in the projection area [compare Fig. 1 (a)], which 
can be seen in the band structure of the (99, 97)-tube very nicely. With increasing sampling, the folded helical vector 
cuts the Brillouin zone more often and hence also closer to the K points. Since the energetically lowest transitions are 
those close to the AT-Point, it is clear that the transition energy becomes smaller for larger tubes. This is a well-known 
result^ and can be seen in our model very easily. 


Systematics of the Z3-mode frequencies for non-coprime CNTs 

We now want to extend the geometrical model presented in the main manuscript to non-coprime tubes, which have 
more than one subband. The situation is shown in Fig. 2 for a (15, 3)-tube, which has three subbands (m = —1,0, -1-1). 
Due to the scattering between equivalent minima and the point symmetry at the F point, the scattering process can 
always be represented by a scattering across the F point at (0,0). Let us now assume that an electron is scattered 
between a minimum on the m = -fl band to a minimum on the m = —1 band, i.e., the electron is scattered between 
—fcmin and fcniin. The resonant phonon wave vector is given by q^es = 2 fc^in. We can now draw the already introduced 
triangles with the same edge points as in the case for coprime tubes. Hence, excactly the same systematics as before 
apply. Furthermore, our geometrical approach instrincally assigns the phonons from the correct band, i.e., with the 
correct quantum number. In our example the electrons gets scattered between the m = 4-1 and m = — 1 band, which 
results in a change of the helical quantum number by Am = —2. Therefore, the phonons must stem from the m = —2 
band, which is equivalent to the m = 4-1 band. Since the rh bands are equidistantly spaced in the reciprocal lattice, 
this result is naturally produced in our model and can be seen very nicely. 
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Figure 2. Illustration of the geometrical model for non-coprime tubes. In this example, we used the helical bands of a (15,3)- 
tube. The numbers -1, 0, +1 represent the helical quantum number fa of each subband. Due to the scattering between two 
equivalent minima and the point symmetry at the F point, the resonant phonon wave-vector is always twice the distance F-fcmin. 
Hence, the same systematics as before apply. 




kx (1/A) 


Figure 3. (a) Family-behavior of the D-mode frequencies for Ef 2 transition. The labels refer to the chiral indices of the outmost 
tube of each branch, (b) Reduced scheme of the geometric model at a single K point, kmin and label the positions at which 
the electronic transition occurs and the resonant phonons stem from, respectively. The solid white lines denote the K — M 
high-symmetry directions, (c) Calculated resonant phonon wave-vectors for the Ef 2 transition projected on a single K point at 
(0,0). The positions of the resonant phonon wave vectors evolve along a spiral into the K point for increasing tube diameter. 


A. D-mode families 

The family-behavior for Raman modes in carbon nanotubes is already known for the radial-breathing modes®’^, where 
each branch is given by the relation 

v = mod(ni — 712)8 = ±1 and /i = 2 ni -|- 712 = const. 

Here, we will show that the D mode follows a similar behavior and that the D-mode frequencies arrange in a Kataura- 
plot-like pattern. Figure 3 (a) shows the calculated D-mode frequencies for the E 22 transition of semiconducting tubes 
(taken from Fig. 3 of the main manuscript). As can be easily seen, the frequencies are arranged in branches, which 
are given by above mentioned relation. The labels in Fig. 3 (a) refer to the chiral indices of the outmost tube of each 
branch. Compared to the RBM-families, the positions of the D-mode branches are flipped, i.e., the v = —1 families 
are on the upper side of the E 22 branch and the i' = -1-1 families on the lower side. This behavior can be understood 
if we reduce our geometrical model to just a single K point. 

If the electronic transition occurs at a distance h and an angle ip near a K point, the resonantly enhanced phonons 
stem from the distance 2h and an angle ip from a K' point (see Fig. 2). Since K and K' are equivalent but rotated by 
180° compared to each other, a phonon with distance 2h and angle (p at K' can be represented by a phonon with the 
same distance but angle p -\- a.t K. Thus, we can reduce the electronic transition and the corresponding resonant 
phonons to a single K point, where the electronic transition and the resonantly enhanced phonon stem from oppsite 
sides. Figure 3 (b) shows the schematics of the reduced model. From the Kataura plot, we know that the electronic 
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Figure 4. Sketch of the restoring forces induced by curvature for the RBM. (a) Model of a nanotube with infinite diameter d. 
All carbon atoms move perpendicular to the sigma bonds a. (b) Model of a nanotube with radius r = d/2 in the position of 
equilibrium (open circles) and after displacement (gray circles). All atoms move in radial direction, (c) Zoom to an atom with 
its closest neighbors. The restoring forces are drawn in red. Figure taken from Ref. 13. 


transitions split up into branches with ly = ±1 families. Thus, also the resonant phonons for the U-mode scattering 
process split up accordingly. This can be also seen in Fig. 3(c), where we plot the resonantly enhanced phonons 
for the E 22 transition. As can be seen, the phonons are arranged in branches forming a spiral that runs into the K 
point with increasing tube diameter. Thus, tubes with small diameters can be found in the outer region, meaning a 
higher H-mode frequency. If we now want to understand the family behavior of the D mode, we have to keep in mind 
that we plot the resonantly enhanced phonons, which are on the iF-point side opposite to the electronic transition. 
Thus, the phonon frequencies behave contrary to the families of the electronic transitions or the RBMs. Since the 
TO phonon branch exhibits a trigonal warping at K that is opposite to the electronic trigonal warping, the splitting 
is not as pronounced as for the RBM families. Moreover, the splitting is barely visible for the Ef^ transition since the 
trigonal warping vanishes in vicinity of the K point. 

For the RBM families, it is well-known that the v = —\ family of the E 22 transition gives a much larger Raman signal 
compared to the v = +1 family®. This behaviour is due to a higher oscillator strength of the optical transition and a 
higher electron-phonon coupling for these tubes®. Although the electron-phonon coupling might be different for the 
TO phonon, the oscillator strength should be the same for the RBM and the U-mode scattering process. Thus, a 
similar behaviour regarding the intensities of the U-mode families can be expected. This would mean that the upper 
family of the U-mode E 22 branch is dominant in an experiment and the lower family gives a much weaker in signal. 


Curvature effect on phonon frequencies in single-walled carbon nanotubes 

In order to understand the dependence of Raman-mode frequencies on the CNT diameter, we introduce a simple spring 
model of a nanotube. The use of a spring model is justified by the results from experiments with isotopes^® and 
calculations using force constant models^’®’^^’^^. The structure of a nanotube is simply approximated by a closed ring 
of carbon atoms. This simple model will be sufficient to simulate the inversely-proportional relation between diameter 
and RBM frequency and the proportional relation between diameter and G“-mode frequency (more precisely, the 
G~ mode related to the TO phonon branch for semiconducting nanotubes). Furthermore, this model will give the 
right magnitude of the frequency shift due to curvature. This leads to a general understanding of the dependence of 
frequency on diameter for in-plane and out-of-plane vibrations. 

A nanotube with an infinite diameter can be approximated by a chain of carbon atoms with a constant atomic distance 
a. The introduction of different a’s would lead to a more precise model, but cannot be handled with this approach. 
For this nanotube with infinite diameter, the RBM is an acoustic out-of-plane vibration and therefore represents a 
translation with zero frequency [Fig. 4(a)]. 
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A nanotube with a given diameter is represented by a closed ring of atoms. Figure 4 (b) shows this for a nanotube with 
radius r and angle a between radius r and atomic bond a. The former, strict vibration orthogonal to the sigma bonds 
a [Fig. 4 (a)] now becomes the breathing vibration of the nanotube, with a decreased angle a compared to a tube with 
infinite diameter. Figure 4 (c) shows the zoom to an atom with its closest neighbors in the position of equilibrium 
(open circles) and after displacement (gray circles). In next neighbor approximation, two effective restoring forces 
(FlesRBM) can be assumed for each atom, marked with red arrows for the upper atom. 

In general, the restoring force F is proportional to the displacement Aa of a spring with F = —fAa, f is the spring 
constant. For this model, Aa can be decomposed into a radial part Aorad and tangential part Aotang- This leads to 
the effective restoring force 

UresRBM = FfesRBM.l + FresRBM,2 = —/(2Aai -|- 2Aa2) (3) 

~ 2/(Aa,.ad,l T Autang,! T Aarad,2 T Aatang,2) (4) 

= 2/[2Aarad(-4)] = ■ 4cos(a) • (-Cy) . (5) 


The effective force constant f* can be found by dividing the effective restoring force FfesRBM by the displacement of 
the atom Ar in radial direction. With Ar cos(a) = Aa it follows 


/RBM(a) = = 4/cos(a)2 . 


( 6 ) 


The effective force constant /* changes with a, which is related to a change in diameter. The relation between a and 
diameter d can be found with the cosine formula 


a{d) = arccos = arccos . 

Insertion of Eq. (7) in Eq. (6) gives the dependence of the RBM frequency on diameter: 


(7) 


WRBM(d) = 1/ShM = 2\l^ ■ - . (8) 

V m \ m d 

This equation reflects the well known f /d behavior for the RBM frequency. 

Based on the same considerations and thoughts for Fig. 5 the effective restoring force of the G~ mode (TO phonon 
mode) can be determined to 


FvesG- = /Aa • 4sin(a) • (e^,) . (9) 

With Atsin(Q;) = Aa the effective spring constant f* can be calculated to: 

/a-(«) = ^^]f|^=4/sin(a)2 . (fO) 

Insertion of Eq. (7) in Eq. (10) and using the relation sin[arccos(a;)] = Vl — the dependence of the G“-mode 
frequency on diameter can be found: 

^G- id) = = 2 ^ 1 . \Jl-Q ■ ( 11 ) 

Without curvature, the frequency of the G~ mode must be equal to the frequency of the G~^ mode (LO phonon 
mode), so 2y^//m can be assigned to 1591 cm“^. The distance to closest neighbors is set to a = 0.213nm, which is 
the mean value for the largest and lowest possible distance of carbon atoms in a graphene hexagon. 

Figure 6 shows the calculated diameter dependence (red) for the RBM and Fig. 7 for the G~ mode, together with 
experimental data (blue). The decrease of the RBM frequency for a diameter range between 0.6 - 1.2nm is 49% 
for the calculations and 47% for the experiment. For the G~ mode, the increase of frequency for a diameter range 
between 0.6 - 1.2 nm is 5% for the calculations and 4% for the experiment. The large difference i n abs olute frequencies 
for the RBM between calculation and experiment can be explained by an overestimation of 2y^f /m = 1591 cm“^. It 
is known that the force constant in graphene (no curvature) is weaker for out-of-plane vibrations than for in-plane 
vibrations^. 

Overall, a very good agreement between experiment and theory is obtained. For this model, the main effect of 
curvature on frequency is a change in the angle a, related to a cosine function for the RBM (out-of-plane vibration) 
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Figure 5. Sketch of the restoring forces induced by curvature for the G~ mode, (a) Model of a nanotube with infinite diameter 
d. All carbon atoms move parallel to the sigma bonds a. (b) Model of a nanotube with radius r = d/2 in the position of 
equilibrium (open circles) and after displacement (gray circles). All atoms move tangential to the tubes wall, (c) Zoom to an 
atom with its closest neighbors. The restoring forces are drawn in red. Figure taken from Ref. 13. 



Figure 6. Calculated dependence of RBM-frequency on diameter (red) from Eq. (8) with a = 0.213 nm and \Jf /m = 795.5 cm 
For comparison, fit of experimental data (blue) from^ is shown. Figure taken from Ref. 13. 


and a sine function for the G~ mode (in-plane vibration). This leads to a specific change of the effective spring 
constant and explains the opposite relation between diameter and frequency for RBM and G~ mode and the weaker 
influence of the curvature on frequency shift for the G~ mode in comparison to the RBM. Since for the D mode the 
atoms move in-plane direction (breathing-like vibration of the carbon hexagons), a similar dependence of the U-mode 
frequency on diameter can be assumed as for the G~ mode. Of course, curvature-induced strain, changes in the 
atomic bond lengths and the rehybridization of sp^ to sp^ orbitals^ might affect the frequency as well, but seem not 
to be the main influence for nanotubes bigger than 0.6 nm. 
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Figure 7. Calculated dependence of G“-niode frequency on diameter (red) from Eq. (11) with a = 0.213 nm and \/f /m = 
795.5 cm“^. For comparison, fit of experimental data (blue) from^*^ Is shown. Figure taken from Ref. 13. 
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